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Abstract. We systematically study the first three terms in the asymptotic expansions 
of the moments of the transmission eigenvalues and proper delay times as the number 
of quantum channels n in the leads goes to infinity. The computations are based on the 
assumption that the Landauer-Biittiker scattering matrix for chaotic ballistic cavities can 
be modelled by the circular ensembles of Random Matrix Theory (RMT). The starting 
points are the finite-n formulae that we recently discovered.— Our analysis includes all 
the symmetry classes /3 G {1,2,4}; in addition, it applies to the transmission eigenvalues 
of Andreev billiards, whose symmetry classes were classified by Zirnbauer— and Altland 
and Zirnbauer.— Where applicable, our results are in complete agreement with the semi- 
classical theory of mesoscopic systems developed, by Berkolaiko et al.— and Berkolaiko 
and Kuipers. 10 ' 11 Our approach also applies to the Selberg-like integrals. We calculate 
the first two terms in their asymptotic expansion explicitly. 
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1. Introduction 

A quantum dot is often modelled by a two-dimensional billiard with holes on the bound- 
ary, whose sizes are proportional to the number of quantum channels in the leads. Because 
of the mesoscopic dimensions of the ballistic cavity, quantum mechanical phase coherence 
plays an important role in the dynamics of an electron inside it. Therefore, the electric 
current has an intrinsic stochastic nature, whose fluctuations are of theoretical and exper- 
imental interest. Furthermore, whenever the classical limit of the dynamics is chaotic, it is 
expected that such fluctuations should be characterized by universal features that are well 
described by Random Matrix Theory (RMT).£Alfiilfiia& 

At low temperatures and voltage the scattering inside the cavity is elastic and is described 
by the Landauer-Biittiker scattering matrix 



1) S:~- 



r t' 

1 mxm b mxn 



t r' 



n x n 



where m and n are the number of quantum channels in the left and right leads respectively. 
The sub-blocks r mxm , t nxm and r' nxn and t' mxn are the reflection and transmission matrices 
through the incoming and outgoing lead. Without loss of generality we shall assume that 
m > n. The dimensionless quantum conductance at zero temperature is given by 

(2) Trttt = Tl + ... + T n , 

Ti, . . . ,T n are the eigenvalues of the transmission matrix tP. Since S is unitary, T 1; . . . , T n 
lie in the interval [0, 1]. 

The Wigner- Smith time delay matrix is defined by 

(3) Q:=^hS- 1 — . 

The eigenvalues T\, . . . , r n of Q are called proper delay times, and their average 

(4) Tw :=-TiQ = -(r 1 + --- + T n ) 

n n 

is the Wigner delay time. The number of proper delay times n is the total number of 
quantum channels in the leads. As the name suggests, r\y is a measure of the time that 
the electron spends in the ballistic cavity. The purpose of this paper is to provide a 
comprehensive asymptotic study in the limit as n — > oo of the moments of the densities of 
the transmission eigenvalues and of the delay times in chaotic quantum dots. 

In two pioneering letters Bliimel and Smilansky— ^ discovered that if the classical limit 
of the dynamics of a scattering system is chaotic, then the spectral correlations of the 
scattering matrix are well described by those of matrices in one of the circular Dyson's en- 
sembles: the Circular Orthogonal Ensemble (COE), the Circular Unitary Ensemble (CUE) 
and the Circular Symplectic Ensemble (CSE). If the dynamics is time-reversal invariant 
and K 2 = 1, where K is the time- reversal operator, then the appropriate ensemble is the 
COE; if K 2 = —1 it is the CSE; if the system does not have any symmetry, then the 
appropriate ensemble is the CUE. 

Under the assumption that the quantum dynamics can be modelled by RMT, the trans- 
mission eigenvalues have the joint probability density function (j.p.d.f.) (Refs. 6.8l l28ll3~4l 
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and HI) 

1 n 

(5) p^ s \T u ... } T n ):=-l[T-(l-T J f 2 J] \T k - T/ } 

j=l l<j<k<n 

where f3 G {1, 2, 4}, a = | (m — n + 1) — 1 and 5 G {—1, 0, 1, 2}. Throughout this paper, 
C refers to a normalization constant that may change at each occurrence. The right-hand 
side of fl5]) is the j.p.d.f. of the eigenvalues of matrices in the Jacobi ensembles: 



C 

3=1 l<j<k<n 



(6) PU*!, . . . ,*„) : = ^n^ /2 (b+l)_l ( 1 - ^) ma+iyi n 

where G [0, 1], j = 1, . . . , n, 

(7) a = — (1 + -)— 1 and b = m — n. 



I _ 1 18 

I -E k ■Ej I j 



The parameter 5 is different from zero only for Andreev quantum dots, which are ballistic 
cavities in contact with a superconductor. For such systems the symmetry classes are not 
Dyson's ensembles, but those associated to symmetric spaces, which were classified by 
Zirnbauer, 74 Altland and Zirnbauer,- and Duenez.— ^ These symmetries are parametrized 
by the pair of integers 

(8) 09, 6) G {(1,-1), (2,-1), (4, 2), (2,1)}. 
The moments of the density of the transmission eigenvalues are 

(9) Mj»(k,n) := (Tr [(tf) k ] ) = (T 1 k + ...+ T n fc > , 

where the angle brackets denote the average taken with respect to the j.p.d.f. (EJ). 

The assumption that the statistical fluctuations for the electric current in quantum dots 
are modelled by RMT leads to an explicit formula for the j.p.d.f. of the proper delay times 
too.— We have 



1 n 

(10) Mlu...,ln):=^U^ /2e ~^ /2 II l7*-7i 



P 



c 

j=l l<j<k<n 



where 7j = r- 1 and th is the Heisenberg time. In our setting th = n. The j.p.d.f. (flOj) is 
the density of the Laguerre ensemble, namely 



n 

;n) p L jx 1) ...,x n ):=in-i /2(6+i) " ie ^ /2 n 



c 

3=1 ^-<3<k<n 







for Xj G [0, 00), j = 1, . . . , n. 

From f JT0|) the moments of the density of proper delay times are the negative moments 
of the density of the Laguerre ensembles: 

(12) M^{-k,n) := -L <Tr Q k ) = 1 + ... + % k ) , k < n/3/2 + 1, 

The above average is taken with respect to the j.p.d.f. (Ilip and the parameter b is given 
by 

(13) b = n-l + 2/(3. 

For convenience, we have used the same notation to denote the parameters b in the Jacobi 
and Laguerre ensembles. In the rest of the article, it will be clear from the context to 
which ensembles it refers to. 
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We compute the first three terms of the asymptotic expansions of the averages (jH]) 
and ( Il2p in the limit as ri — > oo. The starting points of our analysis are the finite- 
n formulae that we computed in a previous publication.— Explicit expressions for finite 
moments usually are rather involved and not suitable for an asymptotic analysis. However, 
our results^ have the advantage of allowing an asymptotic study beyond the leading order. 

Due to their applications to physics and multivariate analysis, the leading order asymp- 
totics of the moments of the eigenvalue density for the Laguerre and Jacobi ensembles have 
been studied in detail both within RMT and using semiclassical techniques.-^"— 
The finite-n moments of the transmission eigenvalues for (3 = 2 were computed by No- 
vaes,— and Vivo and Vivo.— Besides the first part of this work, 53 recently another paper 
appeared 4 - where the finite moments of the transmission eigenvalues for /? £ {1,2,4} were 
computed using a different approach. The leading order term of the density of the proper 



delay times was studied in Refs. I62IJ63I and |68|. However, in these articles our assumptions 
of ideal physical settings do not hold. 

In a pioneering paper Johansson— looked at the fluctuations of the linear statistics 
J2j f( x j) as n — )• oo for ensembles whose j.p.d.f. is 

1 - 

p(xi,...,x n ) = —J^exp(-riV(xj)) Y\_ l^k-Xjf. 

j=l l<j<k<n 

The potential V(x) is a polynomial of even degree and f(x) is a suitable test function. 
He proved that ^ • f(xj) converges to a normal random variable with finite variance. In 
the same paper he computed the next to leading order correction for arbitrary values of 0. 
Subsequently, Dumitriu and Edelman— studied the global fluctuations of the spectra of the 
/3-Hermite and /3-Laguerre ensembles, that is ensembles in which (3 can be any positive real 
number. They also computed the next to leading order correction to the Marcenko-Pastur 
law. Recent results on the global fluctuations of the spectra of /3-ensembles include Refs. 



43 and 19 



The statistical properties of quantum transport have been studied with semiclassical 
techniques too. The semiclassical approach is completely independent of a RMT analysis. 
Although the RMT conjecture for closed quantum systems with a chaotic classical limit 
has a long history,— i ^ JU a a iiLSZ semiclassical techniques have been applied to quantum 
transport relatively recently.- - — i22i2&d£diii£2i£i Within this framework, the elements of the 
scattering matrix are approximated by sums over classical trajectories connecting incom- 
ing and outgoing channels in the leads. What is most important is that the outcome of 
semiclassical calculations are predictions for individual energy-averaged chaotic systems. 
This should be distinguished from the RMT approach, which involves calculating an av- 
erage over an appropriate ensemble of different systems. The semiclassical limit h — > 
corresponds to the limit n — > oo in RMT. 

Semiclassical calculations are perturbative in nature. The inverse channel number n^ 1 
is proportional to h and is the small parameter for the theory. Higher order terms are 
constructed by including successively more intricate families of classical trajectories into 
the sums and then performing the appropriate combinatorics. This task was completed to 
all orders for the conductance^ and for the shot noise.— Higher moments of the density 
of the transmission eigenvalues are also known at leading order from semiclassical theory.- 
Berkolaiko and Kuipers^ 1 ^ computed the leading order generating function of the mo- 
ments of the proper delay times and the first two subleading corrections of the generating 
functions of the moments of both the transmission eigenvalues and the proper delay times 
for (3 = 1 and (3 = 2. Recently, Berkolaiko and Kuipers^ and independently Novaes^ 
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announced two distinct combinatorial treatments of the correlations of the scattering tra- 
jectories that imply exact agreement at all orders in n~ l between the semiclassical and 
RMT calculations of the moments of the transmission eigenvalues. At present, however, 
it seems quite difficult to extract explicit expressions for the asymptotic expansions from 
the combinatorial formulae. The consistency of RMT with the semiclassical analysis of 
mesoscopic systems is a major success of these two approaches. Indeed, in closed systems 
we expect semiclassical theories to be consistent with RMT only when they describe the 
local fluctuations of the energy levels; the averages that we compute are affected by the 
global fluctuations of the spectra of random matrices. It is far from being obvious a priori 
that we should expect an exact agreement at all orders. 

Not many semiclassical results are available for Andreev quantum dots. Adagideli and 
Beenakker- and subsequently Kuipers et al.—^- developed a semiclassical theory of the 
energy gap in the density of states of Andreev billiards. Berkolaiko and Kuipers^ computed 
the first two corrections to the density of states. The leading order term of the density 
of states was previously computed using RMT by Melsen et al.—^ using a perturbative 
expansion developed by Brouwer and Beenakker, 2 - which is based on the computation of 
integrals of matrix elements of matrices in the circular ensembles. 

Since the parameter b in the integrals and (1121) depends on the number of channels 
in the leads, we need to specify how the limit is taken. In the moments both m and n 
tend to infinity, but the scaling parameter 

/ s 771 

14 u = — 

n 

remains finite. In other words, we fix a but let b = n(u — 1) grow to infinity. For the 
Laguerre ensemble it will be convenient to introduce a new variable w by defining 

(15) b = n(w-l) + 2/(3-l. 

The moments (I12p are then recovered by setting w = 2. A similar generalization was also 
considered in Ref. 0. 

Equations. (|7D and (I15p contain n-independent terms which arise due to physical reasons 
and are rather unusual from a purely mathematical point of view. Indeed, they have non 
trivial effects on the subleading order terms of the averages (jHJ) and ffT2]) . (See, e.g., 
Remark 13.81 ) For this reason the moments of the proper delay times cannot be extracted 
from the results of Dumitriu and Edelman--- beyond the leading order. 

When (3 > the density (Q is often referred to as the Selberg density and the averages 



C Jo Jo J f = \ 



as Selberg-like integrals (see, e.g., Refs. I35IJ36I and l50l). Their name originates from the fact 



that by setting k = and C = n, the right-hand side of ffTET) becomes Selberg's integral, 
which was introduced for the first time in quantum transport by Savin and Sommers.— 
They computed the average of the shot noise and the Fano factor for arbitrary (3 nonper- 
turbatively. Subsequently, Sommers et al— and Savin et a/. 65 applied this approach to 
calculate moments and cumulants up to the 4-th order of linear and non linear statistics 
of the electric current. Novaes^ used the Selberg-like integrals to study the moments of 
the transmission eigenvalues for — 2 nonperturbatively. Khoruzhenko et al. 41 combined 
such integrals with the theory of symmetric functions to compute higher order moments 
and cumulants of the conductance and shot noise. 
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One may ask what the asymptotic expansion of ffTBT) is if both a, b — > oo. This problem 



was studied in Refs. I24jj42l and |58| when both a and b are proportional to n. The scaling 
that they used was 

(17) a = (v — l)n and b = (u — l)n, 



with u and v positive constants. The leading order term is known when /3 = 2 (Refs. |24] and 



421 ) as well as for general (3 (Ref. |58|). It turns out that at leading order the integral (1T5]) 
is independent of /3. We compute the first subleading correction for (3 G {1,2,4}. We also 



show that the leading order terms are consistent with the results in Refs. l27J and [58. 

The structure of the article is the following: in Sec. El we outline some of the main ideas 
behind the proofs; in Sec. |3]we treat the asymptotics of the moments for (3 = 2 ensembles; 
Sec. @]is devoted to (3 = 1 and (3 = 4 ensembles; in Sec. [5] we study Selberg-like integrals; 
Sec. E] concludes the paper with comments and open problems. 

Acknowledgements 

We would like to express our gratitude to Gregory Berkolaiko, Peter Forrester, Jonathan 
Keating, Jack Kuipers and Dmitry Savin for stimulating and helpful discussions. 

2. Preliminaries 

As remarked in the introduction, in the first part of this work 53 we computed the aver- 
ages ([9]) and ([12]) for any number of open channels. Our results were expressed in terms of 
finite sums involving binomial coefficients and Pochhammer Symbols, defined by the ratio 

(18) (*)<„, - ^ 

This simple fact has two important consequences. Firstly, the moments ([9]) and ffT2l) admit 
the asymptotic expansions in integer powers of 1/n: 



(19a) n-^fon)"^ 7 * 



oo 



(19b) n k ^M^(-k,n) ~ ^^(^)n^. 

p=0 

Here and in the rest of the paper whenever the parameters a and b appear, they will depend 
implicitly on u or w according to and (fT5]) . 

The second important feature of our results^ 3 - is that the finite-n formulae for and ffT!?]) 
contain meromorphic functions that are particularly suitable to an asymptotic analysis. In 
order to understand this point, it is useful to consider an example for (3 = 2. 

2.1. An Example for (3 = 2. Let us define 

(20) AMjf 6 (fc,n) :=M^ b (k,n)-M^ b (k + l,n). 



u n,k,ji 



In Ref. |53| we found that 

<*> A <^igC)(A) 

where 

ra b (a + b + 2n - 2j + k + 1) (a + b + n) ( k -j+i) {a + n - j + 1) {b + n) [k-j+i) 



(22) a, 



n,k,j 



(a + b + 2n- j)( k+2 )(a + b + 2n- j+ l) {k) (n+ l) ( _ j} 
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The moments are given explicitly by 

(23) M« (k, n) = <(!,„) -gig g) (. £ J t/, 
where 

(24) M«f>(l,„)-^±^ 



a,b 
n,j,ii 



Ja > bK ' 7 a + b + 2n 



is a particular case of Aomoto's integral (see, e.g., Ref. |50j, pp. 309-310). 
We also found a formula for the negative moments in the Laguerre ensemble, 

The moments of the eigenvalue densities in ensembles with orthogonal and symplectic 
symmetries have a similar structure: the summands are products of two factors. The 
first one is n-independent and consists of binomial coefficients. This feature will obviously 
persist to all orders in the asymptotic expansions of the right-hand sides of Eqs. (|2"T|) 
and d25J. 

The coefficients 

appear frequently in enumerative combinatorics, where they are known as Narayana num- 
bers. The n-independent factor in the summands of (125]) are obtained from the Narayana 
numbers by extending them to negative integers: 

<-) ^-iC-OCt-T 1 )- k>0 - 

This identity follows from the elementary relation 

The n-dependent part of the factors in the summands of f[2T]) and f[25l) are ratios of 
Gamma functions, whose asymptotic expansion is^ 

(29) ¥§if) ~ ^B-ii'^fl!"*"!")^, x-^oo, 

where ^^(x) are the generalized Bernoulli polynomials. Now, substituting the leading 
order of ( 129]) into, say, the j-th term in the sum ( 121]) leads to a power of n whose exponent 
is obtained by summing or subtracting the subscripts of the Pochhammer symbols in the 
ratio (122]): 

(30) n i+k-j+i+j+k-j+i-k-2-k+j _ n ^ 

Crucially the exponent is independent of j. Surprisingly this phenomenon happens at all 
orders of n for the moments of both the transmission eigenvalues and proper delay times in 
all the symmetry classes. Thus, at least in principle, to compute the asymptotic series (119aj) 
up to a given order p we need to substitute into (12"T]) the first p terms of the expansion 
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Remark 2.1. In addition to the Narayana numbers we will often need the Narayana 
polynomials 



k 

(3D N k (u): 



k\( k , ; 



as well as their generating function (see, e.g., Ref. Il8l ) 

oo 

p(u,s) := J2N k (u)s k 



(32) 



k=l 



1 - s(u + 1) - VI - 2s + s 2 - 2ms - 2ms 2 + u 2 s 2 

2s 



Indeed our results have a distinct combinatorial flavour. It is interesting to note that also 
the semiclassical approach to these problems relies heavily on combinatorics, although of a 
different kind. From the RMT point of view, it is not obvious a priori why combinatorics 
should play such an important role in these calculations. 

2.2. Generating Functions. The coefficients in the series (If 9al) and flI9b[) become in- 
creasingly involved as we go higher in the order of the expansion. Then, it becomes 
convenient to express our results in terms of generating functions. 

Given the sets of coefficients (u) and V^ p (w), k — 0, 1, . . . , at a given order p in 
expansions (I19a[) and flf9bl) . we define the generating functions 

oo 

(33a) rf)(v):=E^W^ 

k=l 

oo 

(33b) V\?\w,s)-.= Y,V { ff p {w)s k . 

k=l 

If these series are convergent in a neighbourhood of the origin B e , then Dp(w : s) and 
7p s) are analytic in B e and define the moments uniquely. 



As already appeared in the remarks about formula (1211) . our results on the moments of 
the transmission eigenvalues are based on computing the differences between consecutive 
moments. It is therefore convenient to know their generating function too, which is defined 
by 

oo 

(34) &T^ S \u,s):=J2^ p 5 \u)s k . 

k=l 

The coefficients of the powers of s in this formula are defined by 

(35) A7J5«(«):=7if («) -7ffiJ(«). 

The moments at a given order p in the asymptotic expansion (If 9 aft are then obtained from 
the formula 



(36) T^Hu) = T^(u)-Y^^/' M 



8 



The first moments 7^ (w) can be computed directly from the asymptotic expansion of 
Aomoto's integral ( 12~4"|) . which gives 



(37) TZ'°\u 



(u + 1)p +1 



2 (5 



(3 12 



An immediate consequence of the definition (l34"j) is the following functional relationship 
with f!33al): 



08) r p ^\u, s) = 7 ^- 1 (at p ^(u, s) - rZ"\u)) . 

3. Unitary Symmetry 

The techniques and tools that we use are broadly similar for all the symmetry classes. 
However, as is usual in RMT, the power [3 in the Vandermonde determinant of the 
j.p.d.f.'s and f lTT]) affects the complexity of the calculations. Although when 5^0 
in (jSJ) the matrix ensembles that model Andreev billiards do not belong to Dyson's sym- 
metry classes, for the sake of simplicity we shall use the usual terminology when referring 
to an ensemble characterized by a given (3, namely unitary, orthogonal and symplectic if 
= 2, (3 = 1 and (3 = 4 respectively. 

3.1. Leading Order. In Sec. 14. 1[ Remark I4.1[ we shall prove that the leading order 
contributions are the same for all symmetry classes. Therefore, we shall compute the 
limits 

lim n~ x Mf ] (k, n) and lim n fc_1 Mf } {-k, n) 

only when (3 = 2, since in this case the finite-n formulae simplify considerably. 

Proposition 3.1. The leading order terms of the series (119aj) and ( 119bj) are independent 
of (3 and 5 and are given by 



w-l 



( 39a ) v ifi H = { V a - 1 s / /, - 



ifk = l 



J \ U Zl 



i - I (u+ 1)^' +1 



09b) C>)^-£i£(0( 

j=\ J i=l v 7 v 

Proof. We first prove (J39aj. We have 

(40) + = ("")(-*-,) =n i-k w -k- j + 0{n - k) n _, Q0< 

(l + n^-i-j) (l + 7i) ( _ w) 

Inserting the right-hand side of this equation into (f25j) gives 

The infinite series on the right-hand side of (j4ip is convergent for \w\ > 1. Indeed, since 
the ratio of successive terms is a rational function of the summation index j, it is a special 
case of a hypergeometric series. We find 

(«) sgCt-T 1 )G-0""" = ""' iFi(M+ii2i "" ,)> 
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where the hypergeometric function 2-^1(0, b; c; z) is defined in Appendix |B| Eq. (IB.4p . If 
k = 1 this sum is a geometric series which converges to l/(w — 1); if w — 2 it is the mean 
of the Wigncr time delay. When k > 1 applying Eulcr's transformation formula (sec. e.g.. 
Ref. Q, Eq. (15.3.3)) 

2 Fi(a, b; c; z) — (1 - z) c ~ a ~ b 2 -Fi(c - a, c - 6; c; z) 
converts the hypergeometric series in ( l4"2"j) into the finite sum 

A;— 1 

= (u> ! 1)2fc -i - M - *; 2; «r x ) 

1 ^ 

- (w _ 1)2 *-i (j. _ i) (j ._ fc+1) (fc) (j ._ fc+1) (2) (fc _ :> ._ 1) (fc - j - 1)!' 

where we used that (— x)^ = l/(x + Finally, the identity 

1 1 A-1\/Jc-1 

(43) (k - l)( j f_fc + i)(A;)(jf_fc + i)(2)(/ s _j_i)(A; — j - 1)! ~ k^l \ j )\j - 1 

yields the limit (j39aj) . 

In order to prove (139b[) we first set a = and b = (u — l)n in (1221) : 

= ((M + j> ~ 2j + fc + l) — j + l)(j)(^) (fc _j + i) 

/ 44 n n,feJ ~ ((« + l)n -j) (k+2 ){{u + l)n — j + l) (fc) (n + l)^-) 

^ ' u 2k-2j+2 

~ n (^+lp + 0(1) ' 

Inserting this expression into (12 ip and relabelling the summation index gives 

u 2 i 



(45) 



\imn- 1 AM\ 2) (k,n) = -y( k )( k 



(« + !) 



2fe+l ■ 



Finally, Eq. (139b[) follows by taking the limit as n — » 00 of the identity 



fc-i 

(46) n- 1 Mf ) b (fc,n)=n- 1 M J ( ^(l,n)-n- 1 ^AMf ) b (j,n), 

i=i 

where M] 6 (l,n) is Aomoto's integral f )24|) . □ 
An immediate consequence of this proposition is the following. 



Corollary 3.2. At leading order, the generating functions (I33a|) and ( I33b|) are independent 
of (3 and 5 and are given by 



(47a) V ( p >(w,s 



09), „ 5 _ w - 1 - s - v 7 ^" l) 2 -2s(w + l) + s 2 



(47b) ^(v)^ , -ll(. + D. 
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Proof. By definition and using formula ( 13 9 aft we have 



(48) w ~ l h^- 1 ^ 1 

= — —r + -^rp((w, s/(w - If) , 
w — 1 w — 1 v ' 

where Nk(u) are the Narayana polynomials (|3T1) and p(u,s) is their generating function. 
Inserting formula (J3"2j) for p(w, s) gives (I47al) . 

We first compute the generating function A7q^' s \u, s). From the explicit expression of 
the moments (139bj) we obtain 

(49) AT {M (u, s) = -±-p (it 2 , s/(u + l) 2 ) , 

u + 1 

where, as in the proof of Eq. (I47aj) . p(u,s) is the generating function ([52]) . Elementary 
manipulations give 



Finally, Eq. (l47bj) follows from (JMD and (EZJ) . □ 
Remark 3.3. Using Selberg's integral Novaes^ derived the formula 

A ' _J 'k-i\ [2j\ (-iy u j+l 



-eC; 1 ) CO 



n x j J\j J j + i iu + \yi+ 1 ' 

Taking the differences of this sum and comparing to ( |45|) leads to the combinatorial identity 

This formula has already received attention in enumerative combinatorics, where it was 
proved using generating functions by Coker. 26 Subsequently, Chen et al.,— who referred 
to it as Coker's second problem, gave a bijective proof. Combining Novaes' expression 
and (]39bp gives an independent proof of (1521) . 

The identity (|52|) is intrinsically connected to averages in the Jacobi ensemble. Therefore, 
it seems natural to ask whether similar identities should appear when integrating over other 
ensembles, in particular the Laguerre ensembles. It is a straightforward consequence of our 
exact results^ (see also Refs. |32| and |37|) that 

(53) Wm n-^M^tk.n) = \Y ( k \( k )w j . 

In the same paper Novaes also computed this limit using Selberg's integral and showed 
that 

(54) limn^lf 1 ^^ V — ( 2j )( k ~ 1 )w^\l + wf- 2 ^ 1 . 

n->oo » J + 1 V. j J \ 2j J 



The equivalence between the sums ( IB^|) and (|541) appeared in Ref. |25| too, where it was 
referred to as Coker's first problem. 



li 



It is very interesting and surprising that these combinatorial identities can be proved 
solely using techniques from RMT. 



>GB)/„ , 1 - s - Vl - 6s + s 2 



Remark 3.4. When w = 2 the generating function (I47al) reduces to 

(55) vf%a) 2 

which was computed in Ref. [1(3 using both periodic orbit theory and RMT. 

Remark 3.5. The leading order generating function (147bl) agrees with the semiclassical 
formula computed by Berkolaiko et al. 9 Their notation is slightly different from ours. Their 
main asymptotic parameter is the total number of quantum channels N = n+m = (u+l)n; 
they also introduce the variable £ = u/(u+l) 2 . From our point of view, it is more convenient 
to use the dimension of the ensemble as asymptotic variable. 

Remark 3.6. Instead of deriving exact formulae for the moments and extracting their 
asymptotics, one could instead study the asymptotic expansion for the mean eigenvalue 
density 

(56) p n (x) := (y^Jjx-Xj) 

\i=i 

The moments are by definition fx k p n (x)dx. For example, if the average is taken with 
respect to the j.p. d.f of the Laguerre ensemble, then the limit 

rr^ v -i f \ y / (x-z/„)(^ + -x) 

(57) hm n p n (x) = 

n->oo 2nx 

is the Marcenko-Pastur law. — The support of fl57|) is v± = (y/w ± l) 2 , where w is defined 
as in (jl5p . This fact was exploited to obtain the leading contribution to the density of the 
delay times.— The positive moments of this density are known too,— from which the 
corresponding negative moments may be obtained from a substitution x — > x~ x in (|57|) . 
They are easily seen to be in agreement with formula (I39al) . 

3.2. Beyond Leading Order. We now compute the next two terms in the asymptotic 
expansion of the moments of the transmission eigenvalues and proper delay times when 
(5 = 2. 

The starting point to compute higher order corrections to formulae (I39al) and (139bj) is 
to expand the n-dependent factor 

( } (1 + -)(-!-,) 

for the moments of the proper delay times and the coefficients (1221) for those of the trans- 
mission eigenvalues. For this purpose we need the first three terms in the expansion 
(see, e.g., Ref. [l, Eq. (6.1.47)): 



12z 2 



2 ^ (3(a + - l) 2 - (a - (3 + 1)) + 0(z" 3 )^) , z -> oo. 



The coefficients of 1/z and 1/z 2 will give the terms of order 0(n _1 ) and 0(n~ 2 ) respectively 
in the asymptotic expansions ( I19aj) and (119bj) . Such terms will be polynomials in u (or w), 
whose coefficients will be the Narayana numbers (I26p or will closely resemble them. 
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In the Wigner-Dyson symmetry classes there is no contribution at the next to leading 
order when j3 — 2. However, in Andreev quantum dots (5 ^ 0)there is always a non-zero 
correction. 

Proposition 3.7. Let (3 = 2. The next to leading order terms in the asymptotic expan- 
sions (I19al) and f!19bl) are 

(60a) vf^w) = 0, 

where T^ Q '^ (u) is independent of (3 and 5 and is given explicitly in ( I39b|) . 

Proof. For the negative Laguerre moments, we set b = (w — l)n and insert the term of 
order 0{n~ k ) in equation fHUj) into f l25j) leading to 

(6i) vfliw) = Y k ll^V-~i l ) ^ly^i^ + ^ + J + 1 ) -Mi+j)) 



j=0 

oo 

-k-j-i 



(61') _ky^f! ,\ (1' - J-'r 1 

2 3=0 



W 



2 ^ \ k \ k 

3=0 



Equation (I60ap follows from shifting the summation index j — > j — 1 in the sum (I61'|) . 
We now prove (160bj) . We first use fl59l) to expand (T2"Tj) in an asymptotic series and obtain 

where 

FkT( u 'i) = ^l 2 + 1 - - nM / 2 + - J + - J + V2). 
Collecting the coefficients of 5 in (162]) gives 



Elementary manipulations show that the second sum in the right-hand side of (16"5|) vanishes, 
while the first simplifies to 



2(ti + i5^^^72j;vri/2i 



( 64 ) 

5 A / fc \ / fc 
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Finally, the statement of the proposition follows from the combinatorial identity (see 
Ref. H Proposition V7) 



(65) (« + l) 2fc+1 % VLj72j) ([j/2l) nJ ^ (j) (/) J + 1 (ti + 

which allows the polynomial in ( |64|) to be expressed in terms of the leading order term 
T^ 5 \u) given by □ 



Remark 3.8. At first it is quite surprising that the next to leading order term ( I60bj) 



is different from zero. The fact that the correction at order 0(l/n) of the density of the 
eigenvalues of (3 = 2 ensembles is usually zero is a general phenomenon due to the symmetry 
of the ensembles. It was proved by Johansson— for ensembles of random matrices with an 
even degree polynomial potential. It also appeared in the work of Dumitriu and Edelman— 
for the /3-Hermite and /3-Laguerre ensembles. Dumitriu and Edelman— showed that for 
these ensembles, the term of order 0(l/n p ) is multiplied by a palindromic polynomial in 
—2/(3 that is independent of n. At the next to leading order this polynomial is zero at 
(3 = 2. For the Jacobi ensembles this property is still unproved, but given that it is a 
direct consequence of the symmetries of the Jack polynomials, one would expect the same 
behaviour. 

Why, then, is formula (16 Ob}) not zero? The answer is that for physical reasons the 
exponent of the factors (1 — Xj) in the j.p.d.f. of the eigenvalues of the Jacobi ensembles 
does not scale with n (see Eq. ©). If we adopt the scaling ffTT|) . then the next to leading 
order term is zero for (3 = 2. Indeed, this is the content of Proposition 15 . 2 1 for the Selberg- 
like integrals. 

Remark 3.9. For higher order terms, it is no longer possible to simplify the expressions 
which would appear in place of Eqs. (16 ip and (16 2 p through elementary manipulation. 
Instead, one requires a more systematic procedure to handle more complex polynomials, 
which become increasingly involved as we go further in the asymptotic expansion. However, 
as we shall see, often the coefficients in such polynomials are products of two binomial 
coefficients similar to those appearing in the first line of Thus, in many instances 



they can be conveniently expressed in terms of Jacobi polynomials (see Definition IB.ip . 
whose properties can be used to obtain manageable formulae. A first example of this 
procedure is given in the proof of Proposition 13.101 

The next to leading order generating functions will be discussed in Sec. HJ where we shall 
derive a formula for general 5 with (3 G {1, 2, 4}. 

Proposition 3.10. For (3 = 2 the third coefficient in the expansion f!19bp of the moments 
of the proper delay times is 

(66) S«M = ''g' ^W. 



where Pj ■ (x) refers to the Jacobi polynomial of degree j and parameters a and (3 defined 



j 1 

in Eq. (IB.ip and 

w + 1 



(67) w 

w 
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Proof. As in previous proofs, we begin by inserting (j59|) into (|58p and extract the component 
of order 0(n~ k ~ l ). Then, by substituting it into (|25p we arrive at 



where 



,(2)/ • W 2 . _ W 



Gftfaj) = - l)(j + l)(3j + 2) - —(k + j)(k + j + l)j(l + j) 

+ + ^ + J + V(k + 3 + 2)(3fc + 3j + 1). 

The demonstration now proceeds with lengthy and systematic, although elementary, 
algebraic manipulations. It is useful to give an example that clarifies the pattern of the 
calculations. Consider the contribution to (168]) coming from the linear term in G k 2 (w,j): 

_ k(l + kf y>(k + j\ fk + j + l\ t| 



4 ^ U+ 1/ V fc + 1 



Using the definition of hypergeometric function flB.4j) and the relation (IB.5bj) between 
hypergeometric functions and Jacobi polynomials, the right-hand side of fl69|) becomes 



(70) -*'* +1 );'* + 2 » m ->-'(,F l(t + 2, t + 3,2, m -))- + + 



Repeating the same procedure for the other powers of w in G k 2 \(w,j) eventually gives 



24 (w - l) k+3 



(71) iw(A; + l)(fc + 2) ■ v ft ' 

Let us subtract the right-hand side of Eq. (1661) from (!7T|) ; then, the statement of the 
proposition is equivalent to the identity 

(72) 3k(pV> 2 \w) + Pl?' \w))-6kP^ ) <^ 



w 



w 



+ (3* - 2)(P fe ( if (w) - P™(w)) - 2P {2 »{w) = 0. 

The proof follows from the orthogonality of the Jacobi polynomials and the following 
general procedure: 

_____ /2 2") ^_ 

(1) write each Jacobi polynomial in (J72J) in terms of the polynomials Pj ' (w) using 
the connection formulae (1B.3j) ; 

(2) apply the three-term recurrence relation (1B.2[) to each Jacobi polynomial until only 
the polynomials Pj: 2 ' 2 \w) and Pj^(w) appear in the formula; 

(3) the resulting expression is of the form 

A k {w)P^ 2 \w) + B k {w)P^{w\ 
for some rational functions A k (w) and Bu(w); 
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(4) a direct computation shows 

A k {w) = B k {w) = 0. 

□ 

Remark 3.11. Note that step (1) of the above proof amounts to writing each Jacobi 
polynomial in a consistent linearly independent basis. By formula (1B.3j) the number of 
non-zero coefficients in this basis is typically quite small and independent of k. This 
property is essential for our procedure to work. Indeed, consider the identities 



J^M + P« W = (* + 3)(* + 4)^(0) + <* + i)F£»>~ 



w 



(2k + 3)(k + 2) 2k + 3 

3)(k + 4)P^ 2 \w) (k + l)P^l\w) 




(2k + 3)(k + 2) 2k + 3 

which are a direct consequence of flB.3j) . When substituted into f[72]) . they immediately 
yield the desired cancellations. 

In some of the proofs of this paper (more specifically in Lemma 13.131 and Propositions 
13.141 14.8[ 14.91 and 15 .2 j) we will use this method to justify the equivalence of two apparently 
different combinations of Jacobi polynomials. This task is most conveniently performed 
using a computer algebra package such as Maple or Mathematica, whereby the rather heavy 
algebra involved in steps (l)-(4) may be executed algorithmically. 

In order to keep the notation as simple as possible, we now introduce a convention that we 
shall often use when proving statements about generating functions. If p(x) is a polynomial 
or a function analytic in a neighbourhood of x = 0, then we shall write ja^Wa?) to indicate 
the Taylor coefficient of x k in p(x). This formalism can be trivially extended to a function 
of several variables. Now, when in a proof we write [s fc ]7p (it, s) or \s kl \ [ 'D^\w 1 s), we 

do not refer to (u) or V^(w), but to the Taylor coefficients of a function, usually 
an algebraic function, whose explicit expression is the statement we intend to prove, like, 
for example, the right-hand sides of Eqs. (I47aj) and f!47bj) in Corollary 13 .21 In other words, 
adopting this slight abuse of notation, Corollary 13.21 is equivalent to 

[s k ]V^ ( W , s) = v[%w) and [s^* \u, s) = T^\u). 
We now have 

Corollary 3.12. The generating function of the moments fl66|) is given by 



(73) V y 2 '(w,s) 



2 



(s 2 - 2(w + l)s + (w - l) 2 ) 



5/2 ' 



Proof. The idea of proof is straightforward. We first show that the Taylor coefficients of 
f[73l) satisfy a particular three term recurrence relation. Then, we demonstrate that the 
moments fl66l) satisfy exactly the same recurrence relation. Finally, we are left to verify 
the initial conditions 

&)( w 



[«WW) 



[W 



which can be easily checked. 
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Differentiating (J75|) once with respect to s leads to the differential equation 

dV [2 \w,s) 



ds 



s(s 2 - 2(w + l)s + (w - l) 2 )- 

+ (5s(s - u> - 1) + 4(w + l)s - 2(w - If - 2s 2 )vf\w, s) = 0. 
Inserting the power series expansion 

oo 

V {2 \w,s) = Y,Pk(w)s k 

k=l 

and equating coefficients of s k gives 

(74) (w - l) 2 (k - 2) Pk (w) -{w+ l){2k - l)vk-x{w) + (k + l)p k -2(w) = 0. 
We now have to prove that the moments fl66|) satisfy the recurrence relation f !74|) . In- 



serting fl66|) into the left-hand side of f J74j) and simplifying leads to 

(k + 2)(k- 2)p£l\w) - (2k - l)kwP£f(w) + (k — l)kP£l\w), 

which is just the three term recurrence relation flB.2j) for Jacobi polynomials with a = (3 = 2 
and n = k — 3. □ 

As we go beyond the next to leading order in the expansion ( 1 19 aft , the corrections to 
the moments of the density of the transmission eigenvalues become increasingly involved. 
Therefore, although we need the explicit expression of such coefficients in the proof, it is 
more convenient to state our results only in terms of generating functions. 

We first need the following lemma. 



Lemma 3.13. The Taylor coefficients of 
(75) f(s;u) 



u 2 s 



(i- s y/ 2 ((u + i) 2 -s(u-i) 2 ) 5/2 

as a function of s are 



(76) 



where 



~ U2 + 1 

77 u := — -. 

u z — 1 

Proof. The idea of the proof is the same as that one of Corollary 13.121 and consists of two 
parts. First we obtain a recurrence relation for the Taylor coefficients of the generating 
function ( I75P ; then we argue that the right-hand side of (j76| satisfies the same recurrence 
equation. 

One easily sees from (175]) that 

(78) 

[s 2 ]f(s;u) 



u 



(w + 1) 5 ' 

,2/o„,2 



u 2 (3u 2 -4m + 3) 

(u + iy • 
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Differentiating (|75p once with respect to s leads to the differential equation 

(79) 2s(l -s)(A- Bs) df ^' ,U ^ + UBs 2 - 3Bs + As - 2A)f(s; u) = 0, 

as 

where 

A = (u + l) 2 and B — (u— l) 2 . 
Substituting the expansion 



f(s;u) = ^2p k (u)s h 



k=l 

into (1791 and equating the coefficients of s k gives 

(80) 2A(A; - l) Pk (u) + (3A-B- 2(A + B)k)p k ^{u) + 2Bkp k _ 2 (u) = 0, 

with initial conditions given by (ITS]) . 

Inserting the right-hand side of ( ITEj) into the recurrence relation ( 1SD|) shows that the 
statement of the lemma is equivalent to an identity among Jacobi polynomials, which can 
be proved using the strategy outlined at the end of Proposition 13.101 and in Remark 13.111 
For this particular case it is more convenient to express each Jacobi polynomial in the basis 
Pj 1A \u). □ 

It turns out that f(s;u) is the generating function A7^ 2 '°^(m, s). This result was first 
obtained using periodic orbit theory by Berkolaiko and Kuipers— — more precisely, they 
computed ^ 2,0 ^(m, s), which is related to A7^ \u, s) by (j551) . The RMT proof of this fact 
is a particular case of the next proposition, where we obtain the same generating function 
for Andreev billiards. 

Proposition 3.14. Let (3 = 2 and 5 > —2. The generating function for the coefficients 
T$ 5 \u) is 

m _ 5 2 us((u + l) 2 - s(u - l) 2 ) - 4sV 
{ } 2 1 ' J 4(1 - S ) 3 / 2 ((n + l) 2 - S ( M - l) 2 ) 5 /2 • 

Proof. Our techniques to compute the moments of the transmission eigenvalues are best 
suited to compute the differences of the moments. Therefore, we shall study 

(82) ATr(u,s)^ SU2 + Pu{s(U - 1)2 - (U+1) :U 



4 v / T^((m + 1) 2 -s(m-1) 2 ) 5/2 4(u + 1) 3 ' 
Expanding the coefficient (|22|) to order 0{n~ 1 ) and inserting it into (12T|) gives 



< 83 > A7fi«)(«) = iEG%-i)(=T^^)- 

where ( u ^j) is a polynomial quadratic in 5 and of 4-th order in u; the definition of 

Fjff\u,j) is reported in Appendix [Dl Eq. (1D.2j) . Thus, the right-hand side of (IHBl is a 
quadratic form in S too. 

In order to prove this proposition we need to show that Eq. (IHBl coincides with the fc-th 
Taylor coefficient of the right-hand side of (l82|) in a neighborhood of s = and for any 
u > 1. The proof is systematic and we shall outline the main steps. 

Take, for example, 

[S°] [u^F^iuJ) = -{k- 3 ){k- 3 + \){2k - 2j + l)/3 
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and consider its contribution to the right-hand side of Eq. 

2k(k - 1) A (k - 1\ (k - 2\ u 2k ~ 2 i kJ^fk-l\fk-l\ u 2k ~ 2j 



r _ AI , J jxj-1J(u + 1) 2 ^ 3^V j J\j-lJ(u+l) 2k ^ 

The ratio of consecutive terms in these two sums are rational functions of the summation 
indices. Thus, they are special cases of hypergeometric functions, which, in turn, can 
be expressed in terms of Jacobi polynomials using the correspondence (lB.5aj) . Therefore, 
Eq. 0H1]) becomes 

u 2 (u 2 - l) k - 2 



(u + l) 2fc + 3 



(2k(k-l)P^\u)/^ + kP^(u)/3 



Repeating this procedure for each monomial in u and 5 in Ff^j, we see that f l83|) 

becomes a sum of different Jacobi polynomials, leading to an explicit (though complicated) 
formula for each coefficient of S. For convenience we report these formulae in Appendix 0, 
Eqs. flTX3|) . flTXijl and (jESJ) . 

Up to now we have proved that AT k 2 ' (w) can be written as a particular combination 
of Jacobi polynomials. The rest of the proof proceeds as follows: a) we shall express 
[s fc ] A7^ ( ' 2 ' <S ' ) (m, s) as a linear combination of Jacobi polynomials too; b) by using the same 
procedure as in the proof of Proposition 13.101 — points (l)-(4) — it follows that 

[s k ]AT 2 {2 > S \u,s)-AT k f\u) = 0. 

We shall only outline how to achieve a), as we have already discussed how to prove b). 
First note that A7^ 2 ' <5 ' ) ('U, s) can be written in terms of A^ 2 ' ^, s): 

AT«"V. .) = .) ( 1 + - + jr^. 

2 v ' 1 v ' \ 4u 4us J 4(w + l) 3 

Thus, the Taylor coefficients of f[52"j) are 

AT 2 (2 ' 5) («, S ) = AT 2 (2 ' 0) (u, s) + 5 2 [s k ] AT 2 ( 



*5) ,,.,„> X 

4m .[ fl W]A7f^(ti,.)). 



Now, from Lemma [3.131 and Eq. (|82p we have /(it; s) = AT^ 2 (u, s); therefore 
[s fc ]6AT 2 (2,0) (w,s)(w+l) fc+5 



w 2 A;(A; + l)(w-l) 



fe— 2 



1 - u a ) I fc-^Cfi) - 4P£;?(«) - 4uP™(Z 



Inserting this expression into the right-hand side of (185]) gives [s fc ] At£ 2 ' S \u, s) in terms of 
Jacobi polynomials. □ 

4. Orthogonal and Symplectic Symmetries 

Our aim in this section is to compute the first two corrections to the leading order terms 
(139bl) and (139aj) when /3 = 1 and (3 = 4. The finite n formulae for the averages ([9]) and (ITS"]) 
that we obtained— for ensembles with both orthogonal and symplectic symmetries have 
a similar structure. Thus, the techniques and calculations used to study their asymptotic 
behaviour are almost identical for both symmetry classes. Although we will state the 
results for j5 = 1 and j3 = 4, most of the proofs will focus on j5 = 1 ensembles, as they 
contain the essential details of the arguments. 
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4.1. Prom Finite-n to Asymptotics — Preliminaries. The finite- n moments of the 
transmission eigenvalues and proper delay times are 

(86a) M^ b (k,n) = M^ b (k,n-l) 

[k/2} k-2j 



j=l i=0 



and 



(86b) Mil' (~ k ' n ) = M C (~ k ' n ~ 1 

n/2-l n-2j 



+ I Lb (-k,n). 



We give the explicit expressions of S^ 2 ' b ^ 2 (k } n) and S^ 2 (—k,n) in Eqs. (1A.3|) and OA-laO ; 



the terms Ij ab (k,n) and I^ b (—k,n) are defined in (IA.4I) and (lA.lbl) respectively. The 

(2) (2) 

quantities M) [k, n — 1) and (— k, n — 1) are the moments for = 2, whose formulae 
are given in Eqs. (j2"5|) and (|25p . 



Remark 4.1. Consider the asymptotic expansion of S^ 2,b ^ 2 (k, (n — l)/2) and Ij ab (k,n). 
Inserting the leading order of (T211 gives 

(87a) 2^f 6/2 (fc,(n-l)/2)~—- --^ + 0(n- 1 ), „ oo 



l + u) 



2k\ u 2j 



k 

( 87b ) J U^^)-E^J(TT^ + °( n " 1 )' n ~>°°- 

The computation of the right-hand sides requires an analysis similar to that one used 
to obtain (1441) for (3 = 2. It follows immediately that S^ 2 ' b ^ 2 (k,n) and Ij ab (k,n) do not 
contribute at leading order. As a consequence the leading order term in the expansion fl 1 9 a,j) 
is independent of /3 and 5. The same is true for the moments of the proper delay times (I19bl) . 
In other words, the contributions made by the orthogonal and symplectic symmetries of 
the ensembles affect only the higher order terms. This is a general property of the density 
of the eigenvalues in random matrix ensembles. 

Higher order coefficients in the right-hand sides of (j 19 aft and (119bl) are computed by 
determining the contributions of higher order terms in the expansions of s:^ b/2 (k 

S b 2 (k, n), Ij ab (k,n) and It b (~ k,n). Thus, for example, in order to compute a given 
contribution in the expansion of fl86a[) . we shall separately study the series 



. X 

(88a) -M^M-inE^^'' 

p=0 



|fc/2j k-2j , v , v oo 

(««») e e m G + 2i ) •< M2 (k. (» - d/2) ~ e 

j=l i=o ^ ' ^ p=0 



u)n v 
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and 

(88c) h a , b {Kn) -Ys^C 

as n — > oo. It follows that the coefficients in the expansion (I19al) can be represented as 
(89) T^\u) = U^\u) - S^\u) + 3$f{u). 

For convenience, we have expressed the parameters a and b of the Jacobi ensemble in terms 
of (3 and 5 (see Eq. ([7])); thus, in fl89j) we have replaced the superscript J aj f, with ((3,6). 
We shall adopt this notation in the rest of the paper. Similarly, for the proper delay times 
we shall use the notation ((3) instead of L b (see Eq. ffT5j) ). The contributions to other 
quantities can be broken down in the same way. Therefore, we shall write 

(90a) AT^ («) = All£f (it) - AS { *f (it) + AZ^ 5) (it) 

(90b) 2>« («,) = U%(w) - S%(w) + Z« ». 

Without loss of generality, we shall refer to U^\u), AU^p \u) and U_l p (w) as um- 

tary contributions; similarly, we shall call S^ p S \u), AS^ p \u) and S^l (w) symplectic 
contributions. 

4.2. Next to Leading Order — Negative Moments in the Laguerre Ensembles. 

Proposition 4.2. Let (3 G {1,2,4} ; the next to leading order term in the expansion (I19ap 
is 

(91) Vfl{w) - 




j=0 




Proof. When (3 = 2 the right-hand side of (l9~li) is zero, which is consistent with Proposi- 
tion [3T1 we shall only discuss the proof for (3 = 1, as it is identical to that for (3 = 4. 

Take the next to leading order term as n — > oo of formulae ( lA.laj) and (lA.lbj) . Then, 
by proceeding as outlined in Sec. 14.11 we arrive at 



2k + 2] - T 



( 92b) n»=E » ; 

3=0 V J 

To obtain the next to leading order symplectic contribution we need to work a little bit 
more. We obtain 

(93) s%m =f±{^ x ) (* + 2 /_ + ; - 
(930 =>;>;f*^:r 1 )(* + ^; +1 )^- 1 




2k - 1 



i=o 
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The inner sum in line (193'j) was evaluated using Chu-Vandermonde's summation (see 
Example Appendix P . Combining Eqs. flU2ij) . fl92b|) and (JM} into fl90bj) yields 



\ j=o v 7 jf=0 v 



The identity 



3=0 v J/ j=0 v 7 

which is a simple consequence of the binomial theorem, takes care of the first series in ([9 
the second sum is a hypergeometric function: 

£ C 1 j ) 2%v ~ k ~ ] ~ l = w ' k '^m + l,k + l, l, ur 1 ) 



3=0 



The explicit representation of the Jacobi polynomials (IB . 1 1) combined with (I95p completes 
the proof. □ 

Corollary 4.3. The generating function of the next to leading order corrections of the 
moments fl9T|) is 



2\(3 J \s 2 - 2s{w + 1) + (w - 1 

^ (w - 1) ^s 2 - 2s{w + 1) + (w - l) 2 

s 2 - 2s(w + l) + (w- l) 2 

Proof. We now show that the moments (l9Ti) are the Taylor coefficients of the function (196|) . 
Consider the identity 

(97) (f) w3 = \ (C 1 - vM 2fc + (1 + v^) 2fc ) • 

For s sufficiently small we can write 

oo k 



2^^\2j) 4^1- S (l-v^) 2 1-^(1 



s(l + 0^) 



(98) -^i^o— ' • v. + 

l-s(tw + l) 1 
~ 2(s 2 (w-l) 2 -2s(w + l) + l) ~ 2' 

A direct inspection of the definition of the Jacobi polynomials (IB. ip gives the generating 
function 

09) Ieeu) ^^^e^- 1 )*^^)** 

fe=l 3=0 fc=l 

Recall that w = (w + l)/{w — 1). Furthermore, the generating function of the Jacobi 
polynomials Pj* P \w) wh en a = /3 = is (see Ref. [l, Eq. (22.9.1)) 

oo 1 

<«*») S^W-^brr-i- 

K = l 
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Combining Eqs. flHH]) and f TTUOl leads to 

k=l j=0 VJ/ 



k\ 4 k 1 1 

w 3 s K 



2y/s 2 (w - l) 2 - 2(w + l)s + l 2' 



Finally, to complete the proof we introduce the scaling s — > (w — 1) 2 s and subtract 
(HE) from dHHD- □ 

Remark 4.4. When if = 2 and = 1 Eq. (19"6"j) reduces to 

Pj 1) (2, S ) 



1 1 , ^ _ 1 — 3s — yl — 6s + 



2(l-6s + s 2 ) 
which was computed from periodic orbit theory.— 

4.3. Next to Leading Order — Transmission Eigenvalues. The following proposi- 
tion extends formula ( 160bl) to the cases = 1 and (3 = 4. In Ref. H, Eq. (9), they obtain 
the next to leading order correction of the eigenvalue density, though no proof was given. 
For the Wigner-Dyson symmetries, 5 = 0; the corresponding result first appeared in Ref. 



2ll and was derived via integration over matrix elements combined with a perturbative ap- 



proach based on diagrammatic methods. We will present a derivation using a completely 
alternative method. 

Proposition 4.5. For e {1,2,4} and 5 > —2 one has 



2k 



( 102 ) ^ 1 / k 2 k x , 

+ ^(n+ir+ 2 y m2j+i ~ § u o + 1 



u 2 ^ 2 



Proof. As previously we shall discuss only the symmetry 0=1. Let us set a = 5 + 1, 
b = (u — l)n in formula (I86al) . 

Following the general strategy outlined in Sec. 14.11 we see that the unitary contribution 
in (JEHU) is 

where 

F [ kf\^j) = (1 + k - j)(k + 5 - j - 1) + uk(l -5)+ u 2 j(5 - j). 
Simplifying, we find 

which is obtained using the same method which led from (I62p to the equation (1641) . provided 
that the we substitute a = 5 + 1 and n — > n — 1 in 



23 



The symplectic contribution at order p — 1 is obtained by replacing (I87ap into (I88bj) . 
Using the Chu-Vandermonde summation of Lemma IC1[ we obtain 



^ )(u) "lr5C'-0G^)(rT^ 



(104) 

l£4/Y2fc\ fk\ 2 \ u 2 i 



E 



u 2 ^ 2 



^j^\\^jJ \jj J H + u) 2k - 

Finally, we must examine the term Ij ab (k,n) in (|86a|) . Straightforward application of 
( j875jl gives 

(105) &Xkf(u) = 2u(u - l) 2 \u + iy 2k - 2 . 

Inserting (I104p . HI 03 j) and (1 1 5 [) into (I90ap shows that equation f II 2 j) is equivalent to 
proving the combinatorial identity 

u(u-l) 2k = (u+l) 2k+2 AS£f\u) 

*(gC)'-"-soc:. 

where AS^ S) (u) := S^fiu) - Sj^^u) is given in f lTTHj) . 

Equation (I106P can be proved simply by comparing the coefficients of the monomials u 2 ^ 
(w 2jf+1 ) in the polynomials on both sides of the equation. It turns out that the right-hand 
side can be reduced to the left-hand side by repeated applications of Pascal's rule 

k\ f k \ _ (k + 1 
j) + [j + l) ~ [j + 1 

□ 

Using the identity f l6"5"j) . we can write the statement of Proposition 14.51 in terms of the 
leading order f[5Tj) . Thus we have, 

Corollary 4.6. The first correction for the differences of moments can be written in terms 
of formula flBTj) as 



From Corollary 14.61 we can compute the generating function of f)107p . 

Corollary 4.7. The next to leading order generating function of the moments of the trans- 
mission eigenvalues is 



rl M (u,s) 

(108) 



2 \ / us 



f3 J y( s -l)((u + l) 2 -s(u-l) 2 ) 
, & ( + + 1 



2/3 I V ^7 V /( M + 1)2 _ s ( u _ 1)2 s-1 
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Proof. The generating function of the factor 7^ffo( — u ) m ^q. f llOTj) is easily obtained 
from (147bj) . Then, the substitution 

2 



It — 1 

s — > s . 

u + 1 



and some algebra lead to 
(109) 




2/3 ^ Sv /(u + l) 2 -s(u 
Finally, formula (IT08]) follows from (JHHD ■ □ 

4.4. Second Corrections. As we go higher in the order of the expansions (11 9 aft and (119bj) 
the formulae for the coefficients become very cumbersome. Thus, the results are better 
expressed only in terms of generating functions. For the sake of simplicity, we shall only 
discuss the symmetry classes with /3 = 1. 

4.4.1. Negative Moments in the Laguerre Orthogonal Ensemble. 

Proposition 4.8. The generating function of the second corrections to the moments of the 
proper delay times is 



(i), > (w + l)s 2 - (2w 2 - 3w + 2)a +(w- l) 2 (w + 1) 

2 ^ W > S) - S ( S 2_ 2 ( w+ l) s + ( w _ 1)2)5/2 



(110) 

y J (w -l)s + l -w 2 

+ S (s 2 -2(w+l)s + {w-l) 2 ) 2 

Proof. The proof is split into two parts: first we compute the third coefficients in the 
asymptotic expansion of formula (I86bj) : then we show that they coincide with the Taylor 
coefficients of the right-hand side of Eq. (lllOp . 

We start from the coefficients of the unitary contribution U^\(w). Inserting n — > n — 1 
and then b = 1 + (w — l)n in (1231) we find 



<1U) - 2 ^?(*) + 7 <^p£M 

The derivation of this formula is almost identical to the discussion of Eqs. (lo^|) and (I70p in 
the proof of Proposition 13.101 

The term involving the asymptotics of Ii, b (—k,n) is the easiest to compute; it is an 
application of the expansion (1591) to fjATTb]) . Its contribut ion is 

< 112) » = >- ir E(^ +1 , 

The coefficient of the symplectic contribution leads to the double sum 

j=i i=o ^ ' ^ ■■ / 
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where 

G${w, i,j) = (k + i + 2j-l)(k + i + 2j)/2 - j(j - 1) - w{f + z(4j + i + l)/2) . 

Now consider, for example, the contribution from the coefficient of w in G^\(w,i, j). We 
have 

E t i k V-\ ') i k + !- 1 " ') ^ + w + 4 + ^ 

j = l i=o \ / \ / 

= EE (* + * : i" *) (* + i - i " >^( i2 + « - ">< 4 ' + « - j ' + w)- 

J=l i=l \ / \ / 

The finite sum in this equation can be computed exactly in terms of binomial coefficients 
using Chu-Vandermonde's summation techniques discussed in Appendix O The remaining 
infinite series is a combination of hypergeometric sums which may be evaluated in terms 
of Jacobi polynomials using the identity (IB.5bj) . Repeating this procedure and inserting 

the contributions fTTTTj) . (fTP21) and S$(iu) into fl90b|) we arrive at 

2(w - l) k+2 V$(w) = -k(w - 1)P£?H - k(k + l)(P^\iv) + wP^\w)) 

+ (w/6)(w - l)-\k + l)(k + 2) {13P£$(w) - P£f (tS)) 
(114) k + 2 



2(™ - l) fc ~ ^ V2j + 1 



- fc W P fc ( °f (w) - (1 - 3w(A; + l))P fc a i'^ 



w) 



We need to show that the Taylor coefficients of the right-hand side of Eq. (IllOp are given 
by DgH in (HUD. 
First notice that 



V k sr(M + 2\ 

A(w - l) 2k+1 \2j + l) 



k A/2fc + 2\ (w-l)s + l-w 2 

\ / (/' s = S- 



(s 2 - 2(w + l)s + (w - l) 2 y 



which follows from elementary manipulations involving the binomial theorem and geometric 
series. Then, a straightforward application of Corollary 13. 121 gives 

, r k] {w + l)s 3 - (2w 2 -3w + 2)s 2 + (w - l)\w + l)s 

1 ' ^ J {s 2 -2{w + l)s + {w-l) 2 ) 5 / 2 

= - (2w2 - ro 3ro+2) x>i» + ( "- 1) > +1) ie»- 

Insert the explicit formula f )66|) for ^(^O i 11 ^ the right-hand side of Eq. (I116P and apply 
the three term recurrence realation (1B.2j) to remove 

Pi-?^). We obtain 
(u; + l)P fc ( !f (w)(£; + 3) (8w + 7wk - 3w 2 - 3)(Jfe + l)P£f (w) 



A(w - l) k + 2 12(w - l) k + 3 
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Combining this expression with f 1 1 X 5 j) leads to 

(117) 1 ^ Jf=0 

(8w + 7wfc - 3w 2 - 3)(Jfe + l)P fc (2 _f(w) 

+ 12(w - l) fc + 3 ' 

The demonstration that coefficients (11141) are equivalent to the right-hand side of Eq. ( I117p 
involves cumbersome algebraic manipulations, but it is otherwise straightforward. It can 
be carried out systematically following the steps (l)-(4) of the procedure outlined at the 
end of the proof of Proposition 13.101 and further discussed in Remark 13.111 □ 

4.4.2. Moments in the Jacobi Orthogonal Ensemble. As for the second correction of the 
moments of the proper delay times the results for the transmission eigenvalues are better 
expressed only in terms of generating functions. For simplicity, we first discuss the moments 
of the transmission for the classical Dyson's ensemble, i.e. when 5 = 0. 

Proposition 4.9. The second correction of the moments of the transmission eigenvalues 
for (3 = 1 and 5 = is 



lis) r 2 {1 ' 0) (u, s ) 



us 



(s 2 (u-l) 2 + 3us-(u + l) 2 ) 



((w+l) 2 -s(w-l) 2 ) 5/2 (l-s) 3 / 2 



Proof. The proof follows the usual strategy: firstly, we compute the coefficient at p = 2 of 
the expansion (jl9aj) : then we show that they coincide with the Taylor coefficients of f lllSp . 
Furthermore, rather than computing flllSp directly we will look at 

(119) A^%, S )=< S2{U - 1)2 + 3US -[1 +1)2 K U 



((u+i) 2 -s(u-i) 2 f /2 vT^~s (u+iy 



According to f )90ap . to begin with we split AT^°\u) into the contributions of AW/ 1 ' " ' 



A:, 2 



u 



ASfc^ i u ) an d ATj^^ (u) . First we note that straightforward calculation analogous to the 

proof that AT^i(u) = in Proposition 13.71 shows that AX^'^iu) = 0; therefore, we need 
not worry about this term. 

The contribution of the symplectic term involving the double summation in (I88bl) can 
be written as 

[k/2\k-2j , , / , 

where 

Fj$(u, i,j) = (j + i 2 + 4ij + 2j 2 - 2ki + k 2 - Akj) + 2ukj + u 2 {j - i 2 - 2j 2 - Aij). 
For example, taking the coefficient of u 2 in F^l(u,i, j) leads to 

[k/2\ k-2j 



:«» E E © (, + *J «™'0' - * - - w 

j — 1 % — 



J=l 1=1 
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The inner sums in this equation are of the form 

They are Chu-Vandermonde's summations discussed in Appendix O and can be evaluated 
explicitly in terms of binomial coefficients. This allows us to write (I12ip as a sum of Jacobi 
polynomials. Finally, including the contributions of the linear and constant terms in u in 
F^2(u,i,j) we arrive at 

' " ' ' ;! -^"(") = k(k - 1)(1 - u 2 )Ptl\u) + (k/2)(l - 2k)Pg$(u) 



(123) (u - l) k ~ 2 u 

+ (k/2)u 2 P^(u) + k 2 uP^(u), 

where u is defined in (1771) 
Now, (u + l) 2k+3 AT k {1 2 0) 

(u) is a polynomial in u of degree 2k + 1; for convenience we 
separate it into the sum of two polynomials containing only monomials of even and odd 
degree respectively. We write 

(124) (u + l)«^ATff\ u ) = AT^>°(u) + A7#° ),e («) 
Let us focus first on the term containing only odd monomials. We have 

(125) AT^'V) = W^\u) " A< 2 0) '»- 
The contribution of the unitary term gives 

(126) A^ 2 0) '» = -u(u 2 - l) k (k + 1) ((k/6)(u 2 - l)Ptl\u) + (k + l)Pi°'°\u)) . 

This formula is obtained by inserting f[5T?j) into fl2"Tj) . by shifting n — > n — 1, and setting 
a = 1 and b — (u — l)n. 

Since the numerator of u is quadratic in u, the monomials of odd degree in 

can be extracted from (I123p by simple inspection; then, Eqs. (11251) and (11261) give 
- r f_ z ,U = "(*/2)(l - 2A: + 2u 2 + 2ku 2 + u 4 )P£?(fi) 

LL\LL j. J 

(127a) + (« 4 - - l)P£$(u) - k(k + 1)(1 - u 2 fP^(u) 



'l/2)(u'-l)(k + l)(2k + l-u')P^ 



(0,0), ~> 



- (w 2 - l) 2 (k + 1) ((fc/6)(« 2 - l)P^\u) + Pi 
The contribution of A7^ 2 '° (w) can be computed in the same way. We obtain 

AT&0) ' e(M) ~(k- l)((k/2)(u+l)P^\u) + u 2 Ptl\u) - (Jfc/3)P£?(«)) 



m 2 (m 2 — l) fc 2 



+ (fc/6)(5 M 2 - l)P^(fi) + (2k + l)(u 2 - l)P^(u) 
(12Tb) + (2*/3)(* - lyP^iu) - 2fc(fc - 1)(1 - u 2 )P^(u) 

-(k+ 1) 2 (1 - u 2 )p£?(«) - + l)u 2 P^(u) 
+ k(k-l)Ptl } (u). 



28 



Inserting Eqs. (I127al) and (1127bl) into (11241) gives the explicit expression of the second 
corrections of the moments. We now need to demonstrate that they coincide with the 
Taylor coefficients of the right-hand side of (I119p . 

By comparing (11191) with A7^ ( ' 2 ' ' ) ('U, s) in Eq. (I82p we write 

AT 2 (1 ' 0) Ks) = 3AT 2 m (u,s) - A{u+1)2 AT 2 {2fi \u,s) 
( 128 ) " , ^ U , 1V , 



u 

U ' us 



From Proposition 13.141 we know the Taylor coefficients of ATjP'^w, s). Thus, we end up 
with the formulae 

M°AT 2 (1 ' V, s) = m{k + _ {2/m +m + 2)p ( 0)0) 



u(u 2 — l) k 

(129a) -mp^m+W3)(k-i)p^\u 

-(k/2)(k + l)(u-iyu 2 Ptl\u) 
for the component including only monomials of odd powers of u and 

da*) ^tty = {2/3){k + ^ + i)p '- !>(s » - < i/2)(fc + i)p -' ,(a) 

+ 2k(k + l)ifj?(fi) - (2/3)(* - 1)«K> (fl) 

for the monomials of even degree. 

Finally, the equivalence of the pair of Eqs. ( 11271) to (11291) is proved systematically through 
steps (l)-(4) discussed in the proof of Proposition 13.101 □ 

Remark 4.10. Formula (11181) coincides with the generating function computed with semi- 
classical techniques by Berkolaiko and Kuipers.— 

The computation of the second correction for Andreev billiards (5 ^ 0) is almost identical 
to the proof of Proposition 14. 9[ but obviously the algebraic manipulations are much more 
involved. We did not carry out the proof in its entirety. The expression of such a generating 
function, however, can be inferred from general methods and by inspecting lower order 
terms. 

Conjecture 4.11. The generating function of the second correction of the moments of the 
transmission eigenvalues for Andreev billiards with f3 = 1 and 5 > —2 is 



r 2 (1 '%, s) = — — - — ^ — ^ — — 

+ 5s 
+ 5s 
-5 



us 



(s 2 (u-l) 2 + 3us-(u + l) 2 ) 35us 



((u + l) 2 - s(u - 1) 2 ) 5/2 (1 - s) 3 / 2 2 ( u + !) 3 ( s - !) 
((u + l) 2 - s(u - l) 2 ) 2 + 2u5((u + l) 2 - s(u - l) 2 ) 

2((m + 1) 2 -s(m-1) 2 ) 5/2 (1-s) 3 / 2 
(u + l) 2 (w 2 - 5u + 1) - s{u - lf(u 2 - Au + 1) 
2{s - l)(u + l)((u + l) 2 - s{u - l) 2 ) 2 
3s 3 u(u - l) 4 



2(s - l)(w + l) 3 ((u + l) 2 - s(u - l) 2 ) 



2 ' 
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5. ASYMPTOTICS OF SELBERG-LIKE INTEGRALS 

In this section we will be interested in asymptotics as n — > oo of the integrals 



mPm = U.~ IK 

(130) CJ ° J ° V - ' 



(9/2(5+1)-!^ _ y/2(«H-l)-l 



3=1 / 3=1 

— Sjl' 3 dxi • ■ • dx n , k — 1,2, 



x 

l<i<fc<n 



with the assumption that G {1, 2,4} and 

(131) a — (v — l)n and b = (u — l)n, u, v > 1. 

The integrals ( 1130 j) have been referred to Selberg-like by several authors.— Indeed, 
the normalization constant 



C n (/3,a,b) 

(132) 



</o ■/ o i=1 

x jQ \xk — Xj\ dxi ■ ■ ■ dx n 

l<j<k<n 

66 



was evaluated for any (3 > by Selberg in 1944,— who obtained the formula 

(133) cua a, 6) = if rW>w,))r(W2[^u,))r(u(,H) W 2) 

}i r(/J/2( a + 6 + l + n + i))r(l+/J/2) 

The integral (I132p has since been named Selberg 's integral. 

In the first part of this work— we gave finite-n formulae for fl 1 3 [) which imply the 
existence of an asymptotic expansion in powers of 1/n 

oo 

(134) -M { £\u,v)~^M i £l{u,v)n- p 1 n^oo. 

n p=0 

We now compute the first two coefficients in this expansion. 

Expressions for the /3-independent leading order term appeared in Refs. |24| and |42| with 
the assumption that /3 = 2; a different formula is also available in the literature,— whose 
proof was valid for any 0. We introduce a different approach, which gives a simpler 
characterization of the leading order term and allows to compute the next to leading order 
coefficient. 

5.1. Leading Order Term. As for the moments of the transmission eigenvalues, we in- 
troduce the differences 

(135a) AM[ P \u,v) := M { i\u,v) - Mfl^v); 

similarly, for the coefficients of their asymptotic series we write 

(135b) AMgM ■■=Mi%,v)-M^ l!P (u,v). 

Furthermore, we shall set 

M^\u,v) = AM^\u,v) = 0. 
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Proposition 5.1. Let (3 G {1,2, 4}, a — (v — l)n and b = (u — l)n. We have 

136 MrJu,v) = 7 -> . J 7^ TTTT^ • 

M + w J^J^ W V -1 / (w + v) 2 ^ 1 

Proof. It is sufficient to consider /3 = 2 because, as for the moments of the transmission 
eigenvalues, a simple calculation discussed in Remark 14.11 shows that the terms in fl86ajl 
that do not come from the unitary contributions contribute at subleading order. 
The first term in (I136p is the limit of Aomoto's integral, namely 

lim n- 1 M^ ) (l,n)= lim " '" 



n-Kx) Ja - 6 ' ' n-i-oo a + b + 2n u + v 
Inserting the leading order term of the asymptotic expansion (I29p into (I22p gives 

(0 + v)n - 2j + k + l) ((u + v - l)n) {k _. +1) (vn -j + l) U) (un) (k _ j+1) 



U a f> . _ _ 

( 137 ) ^ {i u + v ) n -3) {k+ 2)^ u + v ) n -3 + l ) {k) ^ + l){- j) 

V i u k-J+U u + v - i)fc-i+i 

~n v + 1, n-Kx, 

(u + vy K+1 

where a and b depend on u and v through the scaling (11311) . Substituting the right-hand 
side of ( I137P into fT2TI) leads to a formula for the differences: 

k 



l^{k\f k \ v j u k ~ j+1 (u + v- l) k -i +1 

3=0 



(138) AAO^iEGXA) (u + v y M 



□ 



5.2. Next to Leading Order Term. 



Proposition 5.2. Let (3 G {1,2,4}, a = (v — l)n and b = (u — l)n. T/ie nezi to leading 
order coefficient in the expansion (I134p is 

Proof. We discuss only the case when (3 = 1, as when — 2 and /3 = 4 the proof is 
analogous. 

We can write down a relation analogous to that one in (1891) . where now a = (v — l)n 
grows with n instead of being independent of n. For the first subleading correction it reads 

(139) M[%, v) = U J k ?(u, v) - Sjy (u, v) + 2fc>, v). 

Insert the leading order of (1291) into the symplectic contribution in (I86ap ; then, a calculation 
similar to the derivation of (I104p in Proposition 14.51 shows that 



\k/2\ k-2j . . . . 



fc\ / A; \ 1 + f — l) k % j 



(u + v) 2k 



2fc\ /A;\ \ fu(u + v — 1) 



2(-u + t>) 2fc y V^jV ViV /V 11 
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In the same way, starting from formula (IA.4D we arrive at 

v k sr~^ ( 2k\ / u(u + v 



Therefore, it is sufficient to prove that 

k /, \ 2 



where we have introduced the new variables 

u y + 1 

(140) y — ~ (u + v — 1) and y= . 

u y - 1 

Taking the differences, we see that the statement of the proposition is equivalent to 

(i4D mi rM = -g^r® + ffi^P ffta. 

Replace n — >■ n — 1 in f[2TT) and insert the scaling (11311) . Then, use (159]) to expand (I14ip 
up to the next to leading order correction. We obtain 

^ (M ' V) = k g [j) [j - l) 2(u + v) 2k ^ 

(142) x / (fc-j + l)(fc-j-2) _ J(J + 1) + (fc-j + l)(fc-j-2) 

\ (u + v — 1) V u 

6k + 4 - 2k 2 + Akj 

(W + V) 

which may be written in terms of Jacobi polynomials as 

w. + 2fc i \ , 0jl . 



(143) + 

I (6fc + 4-2fc 2 )(y-l) (ltl) A 
+ 2k(u + v) k ^ {V) ) ■ 

By subtracting (I14ip from f!143j) we obtain an expression involving Jacobi polynomials. 
In order to complete the proof, we need to show that it vanishes. As previously, we 
achieve it systematically, following the steps (l)-(4) outlined at the end of the proof of 
Proposition 13.101 □ 

5.3. Leading Order Generating Function. Let us define the generating function 



(144) H(u,v;s) = J^M { k %u,v)s k 

k=0 



We now present a new proof of the following. 

Theorem 5.3 (Novaes— ). The generating function H(u,v; s) obeys the quadratic equation 

us 1 — s n 

145 H = -. r -. -H 2 . 

u + v — (1 + u)s u + v — (l + u)s 
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Proof. By rearranging summation indices as j — > k — j, we can write formula (I136p in 
terms of a Narayana polynomial of degree k (see ( 1HT|) ): 

= j^^- k g y - J I — ; — J 

f / m(m + u — 1) 
-iVi 



Thus, the generating function of AAiJ?\{u,v) can be easily related to the generating func- 
tion of the Narayana polynomials, which we report in (|32i) . More explicitly 

v ( u(u + v — 1) Sf 



(146) if( W ) := £ AM { k %,v)s k = -^—p I 

k=l U + V \ 

The generating functions H and H are related by 

(147) H = (l - s' 1 ) H + M { ff 
It follows immediately that 



(u + vf 



oK u,v) 



u(u + V — 1) Sf \ It 



(148) if = p , . , 

s — l\u + v \ v {u + v) z J u-\-v / 

The generating function p(x,s) satisfies the quadratic equation (see, e.g., Ref. 

P = (x + p)(l + p)a, 

which leads to 

{us + {s- l)H){s + {s- 1)H) = (u + v){s- 1)H + us. 
Rearranging this expression gives (11451) . □ 

Formula (11 45ft is equivalent (in our notation) to Eq. (58) in Ref. |58| obtained by a 
completely different approach. 

5.4. The Limiting Eigenvalue Density of the Jacobi Ensembles. The leading order 
term of the density (|56|) averaged over the Jacobi ensembles first appeared in the multi- 
variate statistics literature^ 1 ^ and later in the context of chaotic quantum transport, first 
when a = and b = O 6 ^ and then only when a = 0.— More recently it has attracted 
attention in the free probability literature^^I with the same scaling as that adopted in 
this section. It was found that 



fiAn\ f \ r -i f \ ( , W( x ~ X -)( X +~ X ) 

(149) poo(x) = hm n p n {x) = [u + v)- 



2ttx{1 — x 
with support 



d™) A±= [J-^—(l--^—]± J-^—(l 



u + v \ u + v J 1/ u + v \ u + v 



The density (I149j) is normalized to one; in Refs. |23| and |27| a different normalization and 
notation is used. We now show that the leading order contributions of the Selberg-like 
integrals Jwf \u,v) are the moments of the limiting density. 

Formula (11491) makes it apparent why it is helpful to consider the differences between 
the moments, rather than the moments themselves: multiplication by 

x k _ x k+i = ^ _ x } x k 
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cancels the poles in the denominator of (11491) . making the resulting integration tractable. 

Our method is similar in style to that used in Ref. [37| for a similar problem involving 
the Laguerre ensembles. 

Theorem 5.4. The difference of the moments of the limiting density Poo{x) are 

rvd-*)p ( X ) dx = i -T( k )( k )'■'"" 



AM^(u,v). 



Proof. Define 



(152) 



The substitution 



A+ 

k 



X- 
U + V 

2% 



x (1 — x)p oc (x)dx 
J x k ~ l (^Jx{\- + X+j - A_A+ - x 2 ^ dx. 



A- + A + A + — A_ 
x = + cos(fl) 



transforms the integral in the right-hand side of fl 152[) into 

T al(u + v) r , /n\\k-l . 2 //i\ m 

I k = ^- / (cti + a 2 cos(9)) sm 2 (6)d6, 

where a± = A+ "^ A ~ and a 2 = X+ ^ X ~ ■ 
Let p and q and z(9) be such that 

(153) ai = p 2 + q 2 , a 2 = 2pq, z(6) = (p + ge*)*" 1 . 

Then, we have 

|^(6»)| 2 ' 

Set 



Using 



we see that 



sm 



(cti + a 2 cos(6)) 
( = e w z{d) and n = e~ ie z{Q). 
2(g) = l-oos(2g) = l a(l _ eSW) 



a 2 



u + v 8ir 
_ a| 
8tt 



r U(l - e 2W )\z(6)\ 2 d6 

J — 7T 



z{6)\ z d6-m / £(6)ri(6)d6 



These integrals can be evaluated by calculating the Fourier modes of (, n and z and 
applying Parseval's identity. The Fourier coefficients are 
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Parseval's identity gives 



(154) 



^| 2k-2j-2 2j ( k \( k 



Solving p and q in terms of a\ and in Eq. f !153|> and using (I150p allows us to write 



u(u + v— 1) / v 2 4uv(u + v — 1) 



Substituting these expressions into (I154p and replacing j — > j — 1 gives Eq. (I15ip . □ 

6. Conclusions and Outlook 

In this article we have computed the first three terms of the asymptotic expansions in 
the limit as n — > oo of the moments of the density of the transmission eigenvalues and 
delay times in chaotic ballistic cavities. The fundamental assumption is that the chaotic 
dynamics allows us to model the scattering matrix (TjQ) with the circular ensembles from 
RMT. Our formulae are available for all symmetry classes (3 G {1,2,4} — with a few 
exceptions for the second corrections to the leading order terms, for which we did not 
perform the calculations for /3 = 4. For the moments of the transmission eigenvalues we 
treat Andreev billiards too. Finally, we studied the asymptotics of the Selberg-like integrals 
as well. 

Our results on the moments of the transmission eigenvalues and proper time delays 
for (3 = 1 and (3 = 2 symmetry classes agree with those computed using semiclassical 
techniques.- - — 1 ^ It would be interesting to know if the results that we obtain for Andreev 
billiards beyond the leading order agree with semiclassics too. 

In recent announcements Berkolaiko and Kuipers 12 and independently Novaes^ show 
that semiclassical computations lead to the same asymptotic expansions of the moments 
of the transmission eigenvalues as RMT. Their results involve combinatorial expressions 
for the correlations of the scattering trajectories. At this stage, however, it is not clear 
how to extract explicit formulae from the combinatorics. Their formulae include unsolved 
combinatorial problems too. It is a challenging project to go even further and to obtain 
the full asymptotic expansions for the moments of the transmission eigenvalues and proper 
delay times using RMT techniques. Nevertheless, it would be interesting pursue this pro- 
gram, since combinatorics would not appear in the calculations and in the final formulae. 
Furthermore, the equivalence of the two approaches could answer unsolved combinatorial 
problems. 

In principle such asymptotic series can be obtained because the asymptotic properties 
discussed in Sec. 12. II hold at all orders in negative powers of n. The finite-n formulae that 
we computed in the first part of this project— contain products of the form 



(155) n 



djZ + c 



3/ 



One could insert the asymptotic series (f2"9|) into fll55p and repeat the systematic procedures 
developed in this article to find further corrections. In practise, however, when we go 
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beyond the second correction it is not clear how to organise the resulting asymptotic terms 
in a way that would lead to results in closed form, as those obtained in this paper. 

We believe that, as for the asymptotic of the ratio ( }2"9"|) . it may be possible to establish 
recursion relations for the terms in the asymptotic series of (I155p . This would lead to 
recursion relations for the full expansions of the moments of the proper delay times and of 
the transmission eigenvalues. 



Appendix A. Exact Results for — 1 and — 4 Matrix Ensembles 

In this appendix we report formulae for finite-n that we derived in the first part of this 
work- 3 and that are the starting points for our asymptotic analysis. 

Equations (I86al) and (186bj) for = 1 were expressed in terms of certain coefficients which 
are defined below. We also give the corresponding results for = 4. The expressions for 
= 2 are presented at equations fl2H and (1231) . Recall that 



T(x + n) 

0*0 (n) = 



r(x) 

is the Pochhammer symbol. 

In the following we assume that the channel number satisfies n > k0/2. This allows for 
some simplifications; furthermore, a small n is not required for an asymptotic analysis. 

A.l. Moments of the Proper Delay Times. We need the following expressions for the 
Laguerre ensemble: 

(AU] ob , , v _ ( 2b + 27i) ( _ fc _ i _ 2j - +1) (2n -i-2j + % 

where <p ] ^ k n is an exponentially decaying term which has no contribution to the asymptotic 
expansions ( I19bj) for any asymptotic order p. Its explicit expression can be found in Ref. 



53l . Sec. 6.3, Eq. (6.29). The exact moments for = 1 are then given by Eq. (186bj) ; for 



= 4 we have 

M^(-k,n) = 2 k ~ 1 Mil(-k,2n) 

j=l i=0 v 7 v 7 

A. 2. Moments of the Transmission Eigenvalues. The formulae we need for the Jacobi 
ensemble are 

(A.3) 

s a, b(k n) = 2 4 i- 3 (2a + 2n-i-2j + 1) (0 (26 + 2n) (fe _ i _ 2j+1) (2a + 2b + 2n) {k ^ 2j+1) 
M ,U (2n-2j + l) ( _ i) (n + l) ( _ i) (a + n + l) ( _ i) (6 + n) (1 _ i) (a + 6 + n) ( i_ i) 
(2a + 2b + An - 4j + l)(2a + 26 + 4n - 2i - 4j + k + 1) 



x 



(2a + 26 + 4n - i - 2j + l) (1+fc) (2a + 2b + 4n - i - 4j + l) (1+fe) ' 
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and 
(A.4) 

I (k n) =A k V ( 2k \ ^ + & + 2n - 4j - 1 + 2AQQ + b + 7i)) (ifc _ J - ) (i(6 + n))^ 
JaA (a + 6 + 2n-2i-lW ) (i(a + n+l))(_ i )(i(l + n))(_ i) 

The exact moments for /3 = 1 are given in fl86ajl : for /3 = 4 we have 
(A.5) «*,») = Mj'(M) "EE L 



Appendix B. Jacobi Polynomials and Hypergeometric functions 

This appendix contains definitions and identities of Jacobi polynomials and Gauss hyper- 
geometric functions that are needed throughout the paper. They can be found in standard 
references, in particular in the books by Abramowitz and Stegun,- Ch. 15, and Szego,— 
Ch. 4. However, because of their extensive use in many proofs, for the convenience of the 
reader we list them here. 

The Jacobi polynomials may be defined explicitly by the formula 

If a > — 1 and > — 1 they are orthogonal on [—1,1] with respect to the measure 

dp = {l-x) a {l + xfdx 
and satisfy the three-term recurrence equation 

(B.2) 2n(n + a + f3) (2n + a + f3 - 2)P^ P) (x) 

= (2n + a + (3- l)((2n + a + (5)(2n + a + (3 - 2)x + a 2 - 0*)p£f(x) 

-2(n + a- l)(n + /3- l)(2n + a + /3)P^f (x), n = 2,3,..., 

with initial conditions 

P K/3) (x) = 1 and P[ a ' P) {x) = \{a + P + 2)x + \{a-l3). 
They have connection coefficients of the form (p e Z) 



(B.3a) P^\x) = Y. C tv,nP 1 



3=0 
V 



(B.3b) PM( X ) = ^(-^pAT 

j=0 



where 



C 



3,P,n 



p\(a + (3 + n + l) ip -j)(a + (3 + 2n - 2j + 1 + p) 



J J (a + /3 + 2n-j + l) [1+p) (a + n + l) { _ j) 
The Gauss hypergeometric function is defined in the unit circle by the series 

0)i(&)j z* 



(B.4) 2 F 1 (a } b;c;z) := £ 



3=0 x > 3 J 



i 
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it can be analytically continued in the rest of the complex plane. For special values of 
a and b this series truncates and becomes a polynomial. In particular, for a, (3 6 Z it is 
related to the Jacobi polynomials by the identities 



(B.5a) 2 Fi(-n, -n + /?; a + 1; x) = ^ /J - - ^ ' (x - l) n ^P n _ .. 
(B.5b) 2 F 1 (n,n + /3;a + l;x- 1 )= ( ^- ) M^^Hn-u^ 



r(a + l)r(n ^gj-l) _ ^ p(Aa) /x + 1 
+ 1) 1 J 

/3+n 



r(n + a-/3 + l) v ' n ~ p \x-l 
x \ n T(a + l)T(n — a) m + 1 



x - i; r(n) n ~ Q_i - 1 



Appendix C. Generalizations of Chu-Vandermonde's summations 
Many proofs in Sec. H] require computing sums of the form 



(c , b) tc-irOC^-r 1 



where p is an integer. These may be interpreted as generalisations of the classical Chu- 
Vandermonde convolution identities 

(C.2a) r V S ^ /r + S 



i=— m 



(C.2b) J] 



i=—m 



m + ij\l — ij \m + I J 
m + i\ fl — i\ fm + I + 1 



r / \ s J \r + s + 1 



We now present two lemmas which show how the asymptotics of the double sums in 
Eqs. (I86al) and (186bj) are related to sums (IC.la|) and flC.lbj) respectively. We then give two 
examples describing the basic strategy used to evaluate (IC.laj) and (1C.lb|) . 

Lemma C.l. Set 9 = [k/2\ . For any finite set of coefficients {Cjj} we have the identity 

(c.3> Ef(:)G; 2 >^|EG! i )G:>-> 



Proof. Let Sk denote the left-hand side of Eq. (10.3[) . Shifting the summation index in the 
inner sum as i — > k — i — j gives 

(C.4) 3 



i=l j=l 
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where in the last passage we interchanged the order of summation. By decomposing the 
outer sum according to whether 9 — i < or 9 — i > Eq. fl0.4j) becomes 

Now note that the upper limits in both inner sums can be replaced by v. in the first one 
because if j = i + 1, then 

A) -(')-* 

in the second inner sum i > 9 + 1, thus 29 — i + 1 < i. Since j lies in the range 

29 - i + 1 < j < i, 

then 

j + i>29 + l>k and ( . . J = 0, 



+ J, 

yielding no contribution to the sum. Finally, relabelling — > gives the right-hand 
side of Eq. (1(13]) . □ 

Lemma C.2. Suppose that {Cjj} is any set of coefficients such that the double series 

j=l i=o ^ ' ^ 7 

z's absolutely convergent. Then 

OO j / 

-w re -|- 

^ 1 fc- 1 a *; - 1 

i=i i=i \ / v 

Proof. First, shift the summation index i — ?• i — j in (10.50 ; then, (1C.6P is obtained by 
interchanging the order of summation and relabelling the indices. □ 

Lemmas IC.ll and 10.21 are important because the inner sums on the right-hand side of 
(1C.3|) and in (IC.61) can be explicitly computed whenever Cij is a polynomial in i and j. 
This technique is better illustrated with two examples. 

Example C.3. Set p = in (ICTbj) . From (102bl we have 



(c.6) e=EE{ i+ f:!" 1 )f i+ f!!" 1 W j . 



E 

which implies 

j 

E 



re + j - i - 1\ /re + j + i - 1\ /2re + 2j - 1 



re - 1 / \ re - 1 7 V 2Ar — 1 

t=-j 



re + j - z - 1\ A + j + i - A l({2k + 2j-l\ /re + j-T 2 ' 



k-l J V re - 1 7 2 I V 2re - 1 7 V re - 1 

j = l \ / \ / \ \ / \ 

This formula was used to obtain (among others) Eq. (|93|) in Proposition 14.21 
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Example C.4. Now take p = 2 in the sum (IC.lal) . We have 



% 3 G - i) G + i) 1 G - i) G + i ~ V 3 % 3 G - i) G + ')' 

~ k ' % j G - i - l) G + i - l) + kj % 3 G - i - 1 

- j2 S G - G +0 + ^' s C - G 1 7_ i) 



k 

J + i 



l=- 3 t=-j 



Evaluating each sum using the identity (!C2aP and simplifying give 

3 

E 



k \f k \ -2 _ k (Ik -1 



To obtain this formula we used the symmetry of the summands in OC.lap under the sub- 
stitution % — > —i. 

Completely analogous considerations allow the computation of these sums for higher 
values of p. 



Appendix D. An Explicit Example 

Many proofs in this article are based on a systematic approach outlined in Proposi- 
tion l3.10l and further discussed in Remark l3.11[ The algebra involved is rather cumbersome, 
but can be easily done with a symbolic algebra computer package like those contained in 
Maple or Mathematica. In order to give the reader a flavour of what such formulae look 
like, in this appendix we report the explicit expressions that occur in the computation of 
the second correction of the differences of the moments of the transmission eigenvalues 

In the proof of Proposition 13.141 we showed that 

The coefficient Fj: 2 ,f\u,j) is polynomial of 2-nd order in 5 and of 4-th in u. It is given 
explicitly by 



(2,(5)/ - . 5„ 5 2 



(D.2) F k,i = Aj + -Bj + —C, 



j ■ 
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where 



2B j 



(u-1) 



Aj = (j — k)(J — 1 — k) (3f - 6jk — j + k + 3k 2 — l) /6 
+ (j-l-k) (2j-l-2k) (j-k)u/3 

- j (j - 1 - k) (-jk + k + 1 + f - j) u 2 

- j (2j - 1) (j - 1) u 3 /3 + j (j — 1) (3j 2 - 5j + 1) u 4 /Q, 

(j - k) (j - 1 - fc) (2j -l-2k) 

-j(2j - 1) (j - 1> 3 + (1 + 2j) (j - fc) (j - 1 - k) u 
-j(j-l)(2j-2k-3)u 2 



and 



Cj = 1/2 {j - k) (j - 1 - k) + u (k + 1) (j - 1 - fc) 

+ M 2 ^ (A ._ j + !) + ^ + X ) /2 ) -j( k+ l)u 3 +j (j - l)n 4 /2. 



One of the main steps in the demonstration of Proposition 13. 141 was to turn formula fID.lj) 
into a linear combination of Jacobi polynomials. Below we give the final expression. 

Define the transform of a sequence £ = (£z)^ (which may also depend on parameters u 
and k) by the formula 

1 ^ (k\ ( k \ u 2k ~ 2j 

The second correction A7^ 2 ' (u) can be written in terms of the sequences A = (A;)^ , 
= and <7 = (Ci)£o as 

A7g 5) (u) = A/" [A] + ^[5] + ^[C]. 



By the definition of Jacobi polynomials (IB.lj) and minor algebraic manipulations we arrive 
at 

(D.3) (u + l) k+6 Af[A] = k(l -5k + 3k 2 )P { ^(u)(u - l) k - 2 (u 2 + l)u 2 /6 

+ k(k - 1)(2 - 6k)P^l\u)(u - l) fc ~V/6 + (k - l) 2 kP^ll\u)(u - l) fc -V(w + l) 2 /2 
- 2k(k - l)P^\u)(u - l) k ~ 2 u 3 {u + l)/3 + kP^\u)(u - l) k - l u 2 {u + l) 2 
- 2k 2 P^(u)(u - l) k - 2 u 3 (u + l)/3 + 2k(k - l)p£f(u)(u - l) fc ~V(l - u(3k - l)/2). 

The coefficient of | is 

(D.4) (u + l) k+6 N[B] = -k(k - l)p£$(u)(u - l)*- 1 ^ + l)u 2 

+ kP£$(u)(u - l) k -\u + l)u 2 + k(k - l)P fc ( i' 2) (^)(n - 1)*"V = 0. 
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This identity is obtained using formula (IB.3P and the three-term recurrence relation (IB . 2() . 
The coefficient of 4- is given by 

(D.5) (u + l) k+b M[C] = kP^(u)(u - l) fe "V(l + u 2 )/2 



(1 + k)uP^{u){u - l) k -\u + l)/2 + ku 2 P^{u){u - l) k -\u + 1) 



- (fc + l)^^^)^ - + 1) + (* + l^P^'iu)^ - l) k -\u + l)/2 



ifc-i 



\fc-i/ 
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